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SUMMARY 


A  potential-flow  panel  method  has  been  modified  to  calculate  the  effects 
of  a  rotor  wake  on  the  time-averaged  surface  pressure  and  velocity  distribu¬ 
tions  on  a  helicopter  fuselage.  The  rotor-induced  velocities  are  calculated 
by  using  a  vortex-tube  wake  model.  The  calculated  pressure  distributions  are 
found  to  compare  well  with  experimental  data  obtained  from  tests  of  a  wind- 
tunnel  model. 


INTRODUCTION 

The  fuselage  of  a  helicopter  operates  in  a  flow  field  which  is  composed 
of  a  combination  of  its  forward  flight  velocity  and  the  downwash  induced  by  the 
rotor  system.  This  combination  can  range  from  only  the  rotor-induced  field  at 
hover  to  very  little  rotor  downwash  effect  at  high  speeds.  Analysis  of  this 
flow  field  at  any  condition  from  hover  to  high-speed  flight  is  important  in 
designing  the  fuselage  and  its  various  components  (such  as  wings  and  tails) . 

Potential-flow  panel  methods  (ref.  1),  which  use  source  and/or  doublet 
panels  to  model  rather  complex  vehicles,  have  been  developed  to  analyze  the 
flow  field  of  fixed-wing  aircraft.  Much  work  (refs.  2  to  12)  has  been  done  in 
the  helicopter  community  to  analyze  the  induced  flow  of  a  rotor  system.  These 
methods  range  from  simple  actuator  disk  theory  to  unsteady  lifting-surface 
methods. 

Even  though  the  dynamic  nature  of  the  rotor  wake  is  important  in  analyzing 
certain  problems  such  as  vibration,  many  design  details  can  be  analyzed  by 
using  a  time-averaged  downwash  field.  Therefore  a  computer  code  has  been 
developed  to  combine  a  vortex-tube  rotor-wake  theory  with  an  incompressible , 
potential-flow  panel  method.  This  method  calculates  both  on-body  and  off-body 
velocities,  fuselage  surface  pressures,  and  total  loads. 

This  paper  includes  a  discussion  of  the  basic  theories,  which  are  incor¬ 
porated  in  the  present  computer  code,  with  a  general  description  of  the  com¬ 
puter  code  and  the  experimental  test,  which  was  used  as  a  basis  for  validation. 
Experimental  fuselage  surface  pressure  data  were  obtained  from  a  helicopter 
wind-tunnel  model  at  an  advance  ratio  of  0.05.  Comparisons  of  experimental 
data  and  the  analytical  method  are  presented  for  thrust  coefficients  of  0.0034, 
0.0050,  0.0066,  and  0.0082. 


SYMBOLS 

Units  used  for  physical  quantities  defined  in  this  paper  are  given  in  the 
International  System  of  Units  (SI) .  The  positive  senses  of  parameters  are 
shown  in  figure  1. 
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a  Radial  polar  coordinate/R  or  Hake  radius/R  or  control  point 
distribution  (see  fig.  2) 

rotor  radius  or  vortex-ring  radius,  m 
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rotor  side  force 
=  Body  camber/R 

angle  of  attack,  deg  (see  fig.  1(a)) 
angle  of  sideslip,  deg 

rotor-shaft  tilt  angle,  deg  (see  fig.  1(b)) 
vortex  strength 

rotor  inflow  ratio,  (v,,,  sin  a  -  Vj) j VT 
rotor  advance  ratio,  VOT  cos  ct 
air  density,  kg/m3 
source  strength 

-  (d2  '  dl)/(d2  +  dl) 

azimuthal  polar  coordinate,  deg,  or  velocity  potential 


Accession  Jojr 
wris  ojmaI 

DOC  TA1 

Unannounced 

Justification 


®y_ 


f 


Distribution/ 

J&ZV.  Ibi  11  o  j 
j  A  /all  ejjd/or  j 


Dlst 

A 


special 


3 


$  total -velocity  potential 

X  wake  skew  angle  measured  from  vertical  axis ,  deg 

\\>  blade  azimuth  position,  deg,  or  wake-area  contraction  ratio,  A  /A 

(see  fig.  1(b)) 

¥  stream  function 

fl  rotor  rotational  speed,  rad/sec  (see  fig.  1(b)) 

Subscripts : 
f  fuselage 

i  local  condition 

n  normal 

o  at  rotor  disk 

p  perturbation  due  to  fuselage 

R  rotor 
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DISCUSSION  OF  THEORY 

The  numerical  analysis  method  presented  in  this  paper  is  based  upon  the 
panel  method  of  reference  1  and  an  extension  of  the  vortex-tube  rotor-wake 
model  derived  from  references  3  and  4. 


Panel  Method  Theory 

The  panel  method  is  derived  from  the  conditions  that  the  velocity  poten¬ 
tial  $  must  satisfy  Laplace's  equation 

V2$  =  0  (1) 

as  well  as  two  boundary  conditions.  The  first  boundary  condition  is  the 
requirement  of  tangential  flow  at  the  impermeable  body  surface  and  can  be 
written  as 


The  second  boundary  condition  is  that  $  must  approach  the  nondisturbed  free- 
stream  potential  at  infinity.  The  velocity  potential  $  is  comprised  of  the 
free-stream  potential  4JO0»  a  nonuniform  onset  or  rotor-wake  potential  <f>w,  and 
a  disturbance  potential  due  to  the  fuselage  $p.  Therefore  $  can  be 
expressed  as 

*  =  <*>«,  +  K  +  4>p  .  <3> 

The  gradient  of  the  uniform  onset  potential  is  given  in  reference  1  as 

V<Poo  -  <4> 

and  the  gradient  of  the  nonuniform  onset  potential  is  defined  such  that 

V<J>  =  -A?(P)  (5) 

w 

where  Av  is  dependent  upon  the  position  of  point  P.  The  disturbance  poten¬ 
tial  is  given  as  a  source  density  distribution  over  the  body  surface  such  that 
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where  P  is  a  point  in  space  and  g  is  a  point  on  the  surface  which  is  a  dis¬ 
tance  r  from  point  P.  (See  fig.  2.)  The  function  a(q)  is  the  source 
density  distribution. 

Since  equations  (4)  to  (6)  satisfy  Laplace's  equation  and  the  infinity 
boundary  condition,  the  source  function  a  can  be  calculated  by  using  the 
boundary  condition  of  tangential  flow  at  the  body  surface  (eq.  (2)).  For 
point  P  on  the  body  surface,  the  normal  derivative  of  equation  (6)  is  shown 
in  reference  13  to  be 


3<Dn(P) 


=  —2710  (P)  + 


ff  L.  CJLJ 

JJS  3n  [r(P,q)J 


O  (q)  ds 


The  normal  derivatives  of  <t>  and  d>.  are 

oo  W 


=  -n  •  V 


34> 

— ~  =  -n  •  Av(P)  (9) 

9n 

respectively.  Thus,  combining  equations  (7) ,  (8) ,  and  (9)  gives 

2m,<p)  -  JL  h  [tiftst) 0(q)  as  ■  -5  ‘  G». +  15<PII  (10) 

where  the  right  side  of  equation  (10)  is  the  normal  onset  velocity  (uniform 
plus  nonuniform)  and  the  left  side  of  the  equation  is  the  normal  perturbation 
velocity. 

In  panel  methods  equation  (10)  is  typically  solved  by  approximating  a 
body  surface  with  discrete  quadrilateral  panels.  Thus,  equation  (10)  can  be 
expressed  as  a  matrix  equation  of  the  form 

N 

Z  (”i  *  ^ij)aj  =  “5i  *  (^«  +  (i  *  1,  2,  3,  ....  N)  (11) 

j=l  V 


where  N  is  the  total  number  of  panels  and  V^.  is  a  matrix  of  aerodynamic 
influence  coefficients  which  is  a  function  of  geometry  only.  (See  ref.  13  for 
details.)  Therefore  equation  (10)  has  been  replaced  by  a  set  of  linear 
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I 

j 

algebraic  equations.  This  system  of  equations  is  solved  by  standard  numerical 
matrix  techniques  to  obtain  the  source  strengths  Oj . 

The  velocity  vector  at  each  panel  control  point  is  obtained  by  using  the 
source  strengths  and  aerodynamic  influence  coefficients.  These  final  veloci¬ 
ties  are  given  by 

N 

^  *  £  V^Oj  +  (i  =  1,  2,  .  .  .,  N)  (12) 

j=l 


The  panel  pressure  coefficient  cpf  j.  is  defined  as 
_  _  Pi  ~  Pa, 


P/i 


<lco 


(13) 


Since  the  static  pressure  p^  at  the  panel  center  can  be  expressed  as  the 
difference  between  the  total  and  the  dynamic  pressure  at  panel  i 


c  .  pt,i  -  qi 

Cp,X  <ioo 


(14) 


Also,  the  total  pressure  pt  ^  can  be  expressed  as  the  sum  of  the  uniform  onset 
flow  total  pressure  Pt/00  plus  a  local  total-pressure  variation  Apt  i  due  to 
the  rotor  wake.  Therefore 


C  =  Pt'« 
P.i 


^  ~  Pc  +  APt,i 

^oo 


(15a) 


or  since  p  =  q  +  p  , 

t ,  OO  OO  00 
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(15b) 


The  free-stream  static  pressure  p  is  eliminated,  thus  giving 


P*i 
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g» 


(15c) 
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gives 
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and  by  replacing  with  0. SpV^ 


P/i 


<0.5)p|vn>|2  -  (0.5)p|vi|Z  +  Aptfi 
(0.5)  p  |  V J  2  ‘Joo 


(15d) 


which  reduces  to 


(15e) 


The  term  Apt  ^  will  be  discussed  in  the  section  concerned  with  the  rotor- 
wake  modeling.  Body  forces  and  moments  are  calculated  by  summing  the  pressures 
over  the  surface  of  the  body. 


Vortex-Tube  Theory 

The  effect  of  the  rotor  flow  field  is  added  to  the  panel  method  solution 
of  the  body  aerodynamics  as  an  onset  flow  disturbance  Av^  (eqs.  (11)  and  (12)) 
and  as  a  total -pres sure  variation  Apt> (eq.  (15e) )  at  each  panel  control 
point.  The  onset  flow  disturbance  is  determined  from  the  vortex-tube  method, 
and  the  total -pressure  variation  is  determined  from  momentum  theory. 

The  vortex-tube  method  of  wake  representation  which  was  used  is  a  varia¬ 
tion  of  the  method  described  in  references  3  and  4.  In  the  vortex-tube  method, 
the  wake  is  defined  by  a  series  of  concentric  vortex  tubes  each  of  which  con¬ 
sists  of  a  finite  but  large  number  of  vortex  rings.  (See  fig.  3.) 

The  stream  function  at  a  point  relative  to  a  vortex  ring  is  shown  in  ref¬ 
erence  14  to  be 


Y(P)  =  l£(d1  +  d2)[K(T)  -  E(T)l] 


(16) 


where  T  is  the  strength  of  the  vortex  ring,  R  is  its  radius,  d^  and  d2 
are  the  minimum  and  maximum  distances,  respectively,  from  point  P  to  the 
ring  (see  fig.  3),  and  K(T)  and  E(T)  are  complete  elliptic  integrals  of  the 
first  and  second  kinds,  respectively.  The  parameter  T  is  defined  as 


d2  ~  dl 
d2  + 


(17) 
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The  axial  component  of  velocity  of  the  stream  function  ¥  is 


i_  av_ 

Kp  dBp 


(18) 


where  R?  is  the  radial  distance  of  point  P  from  the  vortex-ring  axis.  In 
reference  6,  equation  (18)  is  expressed  as 


V_  =  -I _ (AB  +  CDF) 

n  2TTRp 


(19) 


where 


A  =  K(T)  -  E(T) 

(20) 

d_  -  1 

d„  +  1 

B  =  r  H 

h  r 

(21) 

dl 

d2 

C  =  d^  +  d2 

(22) 

n  =  TE(T) 

(23) 

1  -  T2 

F-1-il 

dr2  +  da2  -  dxd2 

_  (1  +  dr)d!2  -  (1  -  dr)d22 

(24) 

2dr2 

2drd^d2 

r  o 

2tV2 

al  -  Lda  * 

(dr  ~  1)2J 

(25) 

r  o 

9ll/2 

d2  =  K  + 

(dr  + 

(26) 

where  da  is  the  nondimensional  axial  distance  from  the  plane  of  the  ring  to 
point  P  and  dr  is  the  nondimensional  radial  distance  from  the  axis  of  the 
ring  to  point  P. 


The  axial  component  of  velocity  for  a  vortex  tube  can  be  expressed  as 
K 


v»  ■  2  (*m  +  Wi) 

j=l  p'3 


(27) 
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where  K  is  the  total  number  of  vortex  rings.  The  vortex-ring  circulation  T 
is  obtained  by  requiring  that  Vn  be  equal  to  the  average  induced  velocity  at 
the  center  of  the  rotor  disk  with  uniform  loading. 

Reference  4  shows  that  nonuniform  disk  loading  can  be  represented  by  a 
series  of  concentric  vortex  tubes.  Loading  is  either  added  (positive  vorticity) 
or  subtracted  (negative  vorticity)  from  the  tip  vortex- tube  strength  to  vary 
the  loading  across  the  rotor  disk. 

In  tho  present  method,  wake  contraction  was  added  by  varying  the  vortex¬ 
ring  radii  according  to  empirical  relationships  given  in  reference  15.  The 
ratio  of  the  radius  of  the  final  contracted  wake  to  the  rotor  radius  is 
expressed  as 

*00  (  -Ct«58.77\ 

—  =  0.707  +  0.1418^1.0  -  e  )  (28) 

ro 

or  in  terms  of  area 


\ P  = 


(29) 


The  area  ratio  at  a  given  axial  location  zR  is  expressed  as 


=  ip  +  (1.0  -  i|j)eZR  N  (30) 

Ao 


where  N  is  a  slipstream  contraction  factor  based  on  experimental  data  given 
in  reference  15. 

Swirl  is  introduced  as  a  tangential  component  of  induced  velocity  since 
the  lift  vector  of  a  blade  element  is  in  reality  not  perpendicular  to  the  disk 
plane  but  is  perpendicular  to  the  local  resultant  velocity  (ref.  15) .  There¬ 
fore,  the  local  tangential  velocity  Vt  can  be  calculated  by  the  equation 


where  Vj  is  the  average  induced  downwash  at  the  rotor  disk  calculated  from 
momentum  theory. 
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The  pressure  coefficient  used  in  the  panel  method  was  given  as  equa¬ 
tion  (15e)  which  includes  a  total-pressure  term.  The  derivation  of  the  total- 
pressure  term  is  given  in  the  appendix  and  is  expressed  as 


Ap 


t,i  _ 


3oo 


V  +  Atf 
00 

IvJ2 


1.0  + 


(32) 


Equation  (32)  is  evaluated  at  each  panel  control  point. 


COMPUTER  CODE  DESCRIPTION 

Analytical  results  were  obtained  by  using  the  panel  method  of  reference  1 
coupled  wit!  the  vortex-tube  method  as  previously  discussed.  In  the  computer 
code  architecture,  the  wake  calculation  program  was  added  as  a  subprogram  to 
the  panel  method  program.  The  wake  program  is  called  after  all  the  body 
geometry  has  been  input  and  the  control  points  (panel  centroids)  have  been 
calculated .  The  vortex  program  then  calculates  rotor-induced  velocities  and 
total-pressure  variations  at  the  panel  control  points  and  returns  these  values 
to  the  main  program.  Source  distributions  are  calculated  based  upon  the  total 
onset  flow,  and  velocity  distributions  are  calculated  based  upon  these  source 
distributions  and  the  onset  flow.  The  pressure  distributions  are  then  calcu¬ 
lated  based  upon  the  velocity  distribution  and  the  differential  total  pressure 
as  calculated  in  the  vortex- tube  program. 


Modifications  to  the  Panel  Method  Program 

Modifications  to  the  panel  method  program  included  the  addition  of  arrays 
to  store  the  nonuniform  onset  flows  and  total-pressure  values,  addition  of  the 
nonuniform  onset  flow  term  to  equations  (11)  and  (12) ,  and  the  addition  of  the 
total-pressure  variation  term  to  equation  (15e) . 


Vortex-Tube  Program 

The  vortex-tube  program  requires  as  input:  (1)  Coordinates  of  panel  cen¬ 
ters;  (2)  rotor  geometry  including  hub  position,  tip-path-plane  angle  of  attack, 
fuselage  angle  of  attack,  disk  radius,  and  root  cutout;  and  (3)  performance 
parameters  including  Gj.,  VT,  tip  loss  factor,  and  radial  blade  loading. 

Initial  calculations  determine  general  rotor-wake  parameters  which  define 
the  wake  geometry.  These  parameters  include  the  average  induced  velocity,  the 
contraction  ratio,  the  wake  skew  angle,  and  the  static-pressure  differential  at 
the  rotor  disk. 
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The  average  induced  velocity  is  calculated  by  using  the  momentum  theory. 
For  hover 

Vj  -  Vt(o.5  •  C,  •  Aeff)V2  (3 

where  Agff  =  I.oAkij?  -  Rj,2)  is  a  correction  for  the  effective  rotor  area. 
Km  is  a  tip  loss  factor,  and  R  is  the  root  cutout. 


For  forward  flight 


■j(v2  +  X2)3 


Vj  -  0.5  •  VT  •  Cj  ■  Agff 


which  is  solved  in  an  iterative  fashion  since  X  is  a  function  of  Vj.  Equa¬ 
tion  (33)  is  used  as  an  initial  estimate  of  Vj  to  calculate  the  inflow 
ratio  X.  The  inflow  ratio  and  induced  velocity  are  alternately  calculated 
until  Vj  converges  to  a  solution. 

Initial  observations  of  the  skew  angle  measured  from  the  wind-tunnel  test 
(which  is  discussed  later)  indicated  that  the  skew  angle  is  well  approximated 
by  the  expression 


X  *  tan 


1  iw) 


where  \p  is  the  contraction  ratio  and  is  calculated  by  using  equation  (29) . 

The  following  calculations  are  repeated  for  each  panel  control  point. 
After  the  coordinates  of  the  control  point  are  rotated  and  translated  to  the 
rotor  tip-path-plane  (TPP)  coordinate  system,  the  axial  induced  velocity  is 
calculated  by  using  1  to  15  concentric  vortex  tubes  (depending  upon  the  radial 
variation  of  disk  loading) .  Each  tube  consists  of  101  discrete  vortex  rings, 
the  radii  of  which  vary  with  the  expression 


a-[. 


ip  +  (1.0  -  ip)e 


where  zR  is  the  axial  distance  of  the  ring  beneath  the  rotor  and  N  is  the 
empirical  slipstream  contraction  factor.  Therefore,  the  total  axial  induced 
velocity  at  one  control  point  is  the  sum  of  the  contributions  of  101  vortex 
rings  times  the  number  of  vortex  tubes. 

After  the  induced  axial  velocity  is  obtained,  the  tangential  velocity 
(swirl)  is  calculated  by  using  equation  (31) .  The  axial  and  tangential  veloci¬ 
ties  are  then  rotated  into  the  fuselage  Cartesian  axis  system.  The  local 
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total -pressure  variation  is  calculated  by  using  equation  (32) .  These  rotor- 
induced  velocities  and  total-pressure  values  are  then  returned  to  the  panel 
method  program.  These  rotor  computations  add  approximately  5  percent  to  the 
running  time  of  the  panel  method  program. 


EXPERIMENTAL  DATA 

Validation  of  new  analytical  methods  generally  requires  experimental  data 
as  a  basis  for  comparison.  Although  reference  16  presents  surface  pressures 
of  three  helicopter  models  with  a  1 -meter-diameter ,  two-bladed  rotor,  the 
pressure  orifice  locations  were  considered  too  sparse  to  allow  for  accurate 
comparisons  with  the  analytical  method.  Therefore  a  wind-tunnel  investigation 
using  a  helicopter  model  with  more  closely  spaced  pressure  orifices  was  con¬ 
ducted  in  the  Langley  V/STOL  Tunnel.  A  photograph  of  the  model  in  the  test 
section  is  presented  as  figure  4. 


Model  Description 

The  rotor  system  used  in  this  investigation  consisted  of  a  3.15-meter- 
diameter,  four-bladed  rotor.  The  hub  was  fully  articulated  with  flapping  and 
lagging  hinges  coincident  at  4.8-percent  radius.  The  rotor  blades  had 
-8°  twist,  an  untapered  planform,  and  square  tips.  Details  of  the  rotor  system 
may  be  found  in  table  I. 

The  fuselage  shape  for  this  test  is  mathematically  defined.  At  a  given 
fuselage  station  x,  the  cross-section  y-  and  z-coordinates  are  defined  by  the 
local  fuselage  height  H,  width  W,  camber  line  ZO,  and  elliptical  power  N 
as  follows. 


The  parameters  H,  W,  ZG,  and  N  are  obtained  by  applying  the  function 


f(x)  =  C6  +  C7 


Ci  +  c2 


KT 


Vc8 


(37) 


with  separate  sets  of  constants  to  C8  for  each  of  the  four  cross-section 

parameters  H,  w,  ZO,  and  N.  As  shown  in  figure  5,  the  fuselage  is  divided 
into  four  regions  and  the  pylon  into  two  regions  with  a  set  of  constants  for 
each  region.  These  constants  are  given  in  table  II. 

By  using  the  three  parameters  H,  W,  and  N  as  obtained  from  equa¬ 
tion  (37) ,  the  cross  section  can  be  defined  by  the  polar  coordinates  r  and  <t> 
such  that  for  a  given  $> 


r(4>) 


JO. 


(0.5H  +  0.5W) 


N 


5H  sin)N  +  (0.5W  cos) 


1 


1/2 


(38) 
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The  body  Cartesian  coordinates  are  then  obtained  from  the  relationships 
y  =  r  (4>)  sin  <j> 
z  =  r(4>)  cos  4>  +  ZO 


(39) 

(40) 


Test  Conditions 

The  experimental  data  presented  in  this  paper  were  obtained  for  an  advance 
ratio  of  0.05  at  the  test  conditions  as  given  in  table  III.  The  wind-tunnel 
test  section,  which  measures  4.42  meters  by  6.63  meters,  was  configured  with 
the  walls  and  ceiling  removed.  Angle  of  attack  and  test-section  dynamic  pres¬ 
sure  were  corrected  for  wall  effects  by  using  methods  described  in  reference  17. 
All  tests  were  conducted  with  the  moment  reference  center  of  the  model  on  the 
center  line  of  the  test  section.  This  position  corresponds  to  a  rotor  height 
of  1.84  times  the  rotor  radius.  The  pressure  data  were  averaged  for  50  samples 
over  5  seconds  of  sampling  time.  The  data  used  in  this  paper  were  obtained  from 
reference  18. 


PRESENTATION  OF  RESULTS 

The  results  obtained  from  the  computer  code  and  the  experimental  tests  are 
presented  in  the  following  figures: 


Figure 


Isolated  rotor-induced  velocities  in  y  =  0  plane  for  y  =  0.05  ...  6 

Calculated  x-  and  z-components  of  total  velocity  on  left  side  of 

fuselage  for  y  =  0.05 .  7 

Pressure  distribution  along  top  center  line  of  fuselage 

for  y  =  0.05  8 

Lateral  pressure  distributions  at  four  fuselage  stations 

for  y  =  0.05  9 

Separation  and  swirl  effects  on  fuselage  cross-section  pressure 

distribution  .  10 


DISCUSSION  OF  RESULTS 

Results  obtained  from  the  analytical  method  are  presented  in  the  form  of 
rotor-wake  onset  velocity  distributions,  total-surface  velocity  distributions, 
and  surface  pressure  distributions.  These  three  outputs  are  useful  in  analyzing 
the  aerodynamics  of  a  helicopter  in  different  ways. 

Figure  6  presents  the  induced  velocity  vectors  of  the  isolated  rotor  in  a 
vertical  plane  along  the  center  line  of  the  hub  (i.e.,  yR  =  0)  at  four  thrust 
levels.  Evident  in  this  figure  is  the  contraction  and  skew  angle  of  the  rotor 
wake  as  calculated  by  the  rotor-wake  segment  of  the  program.  As  would  be 
expected,  the  skew  angle  decreases  with  increasing  thrust  as  the  momentum 
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r-  f 


inflow  velocity  increases.  The  skew  angle  ranges  from  approximately  40°  at  the 
low  thrust  to  25°  at  the  high  thrust.  Full  wake  contraction  occurs  at  approxi¬ 
mately  50-percent  radius  below  the  rotor  disk. 

The  inflow  variation  at  the  rotor  disk  is  also  evident  with  zero  inflow  at 
the  center  and  maximum  downwash  at  approximately  80-percent  radius.  There  is 
also  a  substantial  difference  in  inflow  levels  for  the  front  and  aft  portions 
of  the  rotor.  These  variations  decreased  downstream  of  the  rotor. 

In  figure  7  the  x-  and  z-components  of  the  total-velocity  vectors  on  the 
left  side  of  the  fuselage  are  plotted.  These  results  are  obtained  from  the 
panel  method  as  rotor-wake  effects  on  the  fuselage.  Wake  boundaries  are  dis¬ 
cernible  on  the  nose  and  tail  cone  as  large  gradients  in  the  vertical  component 
of  the  total  velocity.  These  total  velocities  follow  the  same  general  pattern 
of  the  isolated  rotor-induced  velocities.  Longitudinal  variations  of  vertical 
velocity  are  particularly  evident  on  the  tail  cone.  The  effect  of  the  root 
cutout  is  seen  as  a  stagnation  region  below  and  slightly  behind  the  hub  loca¬ 
tion  at  35-percent-body  length. 

Figures  8  and  9  present  comparisons  of  surface-pressure  coefficients  for 
experimental  data  and  results  obtained  from  the  analytical  method.  Although 
the  pressure  distribution  on  the  entire  fuselage  surface  is  not  shown,  these 
figures  illustrate  the  character  of  the  longitudinal  and  lateral  pressure 
distributions . 

In  figure  8  the  pressure  coefficients  and  geometry  of  the  top  center  line 
of  the  fuselage  are  given.  Since  the  top  center  line  would  be  considered  the 
stagnation  line  for  a  rotor  downwash,  the  large  positive  values  of  Cp  should 
be  proportional  to  the  square  of  the  downwash  velocity.  Very  evident  in  this 
figure  is  the  nonuniformity  of  the  longitudinal  downwash  distribution  as  is 
seen  in  the  velocity  vectors  of  figures  6  and  7. 

Even  though  the  overall  character  of  the  experimental  and  analytical  dis¬ 
tributions  match,  several  discrepancies  are  evident.  For  the  low-thrust  case 
(fig.  8(a)),  the  wake  impingement  location  (point  of  maximum  Cp)  is  not 
coincident  (x/R  ~  0.25)  because  the  analytical  model  assumes  a  constant  skew 
angle  and  the  actual  wake  is  curved. 

Other  discrepancies  in  figure  8(a)  include  the  overall  levels  of  the  down- 
wash  with  the  analytical  method  underpredicting  the  pressure  coefficients  and 
the  seemingly  higher  experimental  downwash  from  the  front  part  of  the  rotor 
disk  as  compared  with  the  aft  portion.  These  discrepancies  are  not  evident  at 
the  three  higher  thrust  levels,  although  the  analytical  method  does  slightly 
underpredict  the  pressures  on  the  aft  portion  of  the  fuselage.  This  underpre¬ 
diction  may  at  least  be  partly  due  to  the  probable  separation  of  the  rotor 
pylon  which  is  not  modeled  in  the  analytical  method. 

Experimentally,  the  forward  location  of  impingement  varied  from  x/R  =  0.25 
for  CT  =  0.0033  to  x/R  =  0.15  for  Op  =  0.0082,  which  corresponds  to  skew 
angles  of  approximately  40°  to  25°,  respectively.  The  aft  impingement  point 
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cannot  be  discerned  from  the  experimental  data  as  it  occurs  slightly  aft  of  the 
last  pressure  orifices. 

Figure  9  gives  the  pressure  distribution,  experimental  orifice  location, 
and  panel  centroid  locations  at  four  fuselage  stations  and  four  thrust  levels. 
These  four  stations  are  x/R  =  0.20,  0.30,  1.34,  and  1.53.  Two  are  ahead  of 
the  rotor  hub  and  two  are  behind.  The  experimental  data,  especially  at  the 
higher  values  of  Op,  in  these  figures  generally  follow  am  M-pattem  with  large 
negative  values  of  Cp  where  the  downwash  accelerates  around  the  cross-section 
comers  and  large  positive  values  of  Cp  at  the  top  (z/R  maximum)  and  bottom 
(z/R  minimum)  stagnation  points.  The  analytical  data  have  this  same 
M-character  but  quantitatively  differ  at  two  regions.  The  analytical  method 
seems  to  underpredict  the  flow  acceleration  around  the  top  corner  (z/R  posi¬ 
tive)  and  to  overpredict  the  stagnation  pressure  on  the  bottom  side.  The 
reason  for  the  first  discrepancy  may  be  the  coarseness  of  the  paneling  at  the 
corners  and  the  fact  that  the  wake  program  generally  underpredicts  the  downwash 
levels.  The  second  discrepancy  is  attributed  to  flow  separating  at  the  lower 
corner  of  the  cross  section  as  illustrated  in  figure  10.  As  mentioned  pre¬ 
viously,  a  separation  model  was  not  included  in  this  panel  method  program. 

Also  evident  in  figures  9  and  10  is  the  effect  of  swirl.  A  crossflow 
on  a  fuselage  section  would  decrease  the  flow  acceleration  at  the  corners  on 
one  side  and  increase  the  acceleration  on  the  opposite  side;  for  example, 

CT  =  0.0082.  At  x/R  =  0.20  (see  fig.  9 (m) ) ,  the  difference  in  the  experi¬ 
mental  pressures  at  the  top  corner  for  the  right  and  left  side  is  4.75  with 
the  left  side  being  higher  (more  negative  Cp) .  The  analytical  results  again 
underpredict  these  corner  pressures.  Behind  the  rotor  at  x/R  =  1.53 
(fig.  9 (p) ) ,  the  pressure  orifice  spacing  was  not  fine  enough  to  indicate  the 
quantitative  difference  of  the  right  and  left  side,  but  the  right-side  pressure 
distribution  indicates  a  higher  acceleration.  Therefore  the  crossflow  right 
to  left  forward  of  the  hub  and  left  to  right  aft  of  the  hub  indicates  a  defi¬ 
nite  swirl  in  the  direction  of  rotor  rotation  as  modeled  by  the  wake  program. 


CONCLUSIONS 

The  surface  pressure  and  velocity  distributions  about  a  helicopter  fuse¬ 
lage  immersed  in  a  rotor  wake  have  been  calculated  by  using  a  potential-flow 
panel  method  having  a  vortex-tube  rotor  model.  The  pressure  distributions  have 
been  compared  with  experimental  results  obtained  from  a  3.14-meter-diameter 
rotor  wind-tunnel  model  with  a  general  helicopter  fuselage.  The  conclusions 
are  summarized  as  follows; 

1.  The  measured  and  predicted  pressure  on  the  fuselage  at  low  rotor-wake 
skew  angles  varies  significantly  both  longitudinally  and  laterally  as  a  result 
of  the  influence  of  the  rotor  wake. 

2.  The  panel-method— vortex-tube  combination,  which  accounts  for  skew 
angle,  contraction,  and  total-pressure  variations,  can  adequately  calculate  the 
character  of  the  time-averaged  pressure  and  velocity  distributions  if  the  tube 
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geometry  (skew  angle  and  contraction)  is  modeled  correctly,  the  total-pressure 
corrections  are  accounted  for,  and  the  swirl  velocities  are  added. 

3.  Increased  accuracy  would  be  obtained  with  the  inclusion  of  a  flaw- 
separation  technique  and  a  curved  rotor-wake  model. 


Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  VA  23665 
April  17,  1980 
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APPENDIX 


CALCULATION  OF  TOTAL  PRESSURE  IN  THE  ROTOR  WAKE 

Momentum  theory  in  hover  considers  the  rotor  as  a  disk  as  shown  in  the 
following  sketch: 


Location  0  is  at  infinity;  location  1  is  immediately  before  the  disk;  loca¬ 
tion  2  is  immediately  following  the  disk;  and  location  3  is  at  the  point  of 
full  contraction.  The  average  induced  velocity  through  the  disk  is  Vj,  and 
the  fully  contracted  wake  velocity  is  w.  The  static  pressures  are  pQ,  p^, 
P2,  and  p3  at  locations  0,  1,  2,  and  3,  respectively. 

Bernoulli's  equation  for  steady,  inviscid,  incompressible  flow  with  negli¬ 
gible  gravitational  force  and  no  work  done  is 

pt  =  p  +  0.5pv2  (Al) 

where  pt  is  the  total  pressure,  p  is  the  static  pressure,  p  is  the  den¬ 
sity,  and  V  is  the  local  velocity.  Applying  equation  (Al)  to  the  four  loca¬ 
tions  and  assuming  that  pt  Q  =  Pt,l  and  Pt,2  =  Pt,3  gives 

P0  =  Pi  +  O.SpVj2  (A2) 

and 

P3  +  0.5pw2  *  p2  +  0.5pVj2  (A3) 

The  thrust  on  the  disk  can  be  expressed  as 

T  *  (p2  -  Pi)A  (A4) 

where  A  is  the  disk  area. 


gives 


Subtracting  equation  (A2)  from  equation  (A3)  and  solving  for  (p2  -  p^J 
P2  "  Pi  "  P3  "  P0  +  0 • 5pw2 
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However  P3  ”  Po  “  °»  the  assumption  is  made  that  the  static  pressure  in 
the  fully  contracted  wake  is  equal  to  the  ambient  static  pressure.  By  making 
this  assumption  and  combining  equations  (A4)  and  (A5) ,  the  rotor  thrust  is 


expressed  as 

T  =  0.5pw2A  (A6) 

From  momentum  theory 

T  =  pAVjW  (A7) 

Equating  equations  (A6)  and  (A7)  and  solving  for  w  gives 

w  =  2 .OVj  (A8) 

Using  equation  (A8)  in  equation  (A5)  gives 

P2  "  Pi  =  2-OpVj2  (A9) 

From  equation  (A3) ,  assuming  that  p2  =  p0 

P2  "  P0  =  0 . 5pw2  -  O.SpVj2  (A10) 

Again,  by  substituting  w  =  2.0V-J.,  equation  (A10)  becomes 


P2  -  P0  =  1-SpVj2 


(All) 


or 


p2  -  p0  =  0.75(p2  -  px)  (Al2) 

Therefore  the  static-pressure  rise  immediately  behind  the  disk  is  three-fourths 
of  the  static-pressure  rise  across  the  disk.  By  applying  the  definition 


c  =  T 
uT  - 


PAVT' 


to  equation  (A4) ,  P2  “  Pi  becomes 


P2  “  Pi  *  CijpVj 

Substituting  equation  (A14)  into  equation  (A12)  gives 
P2  "  Po  "  0.75&PPV.J,2 


(A13) 


(A14) 


(A15) 
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or  by  nondimensionalizing  by  O.SpV^  and  making  the  approximation  y  ■  V00/VT, 
equation  <A15)  becomes 


P2  ~  Pp  =  1-50;. 
0*  SpV^  y2 


(A16) 


which  expresses  the  difference  in  static  pressure  across  the  wake  boundary. 


The  difference  in  total  pressure  inside  the  rotor  wake  at  a  point  i  and 
free  stream  can  be  expressed  as 


APt,l  _  APj  A<U 

<3co 


(A17) 


By  assuming  that  the  static-pressure  variation  at  point  2  in  forward  flight  is 
the  same  as  hover  and  that  the  static-pressure  variation  goes  to  zero  for  the 
fully  contracted  wake,  p  can  be  expressed  as 


where  A^  is  the  wake  cross-section  area  at  point  i. 

The  difference  in  dynamic  pressure  can  be  expressed  as 


Aq±  _  |v.  +  avJ2 
"  IvJ2 


(A19) 


where  Av^  is  the  resultant  of  the  axial  and  tangential  wake  velocities. 
Therefore 


APt,i 

<Ico 


|V,  +  Avil 

I  v_| 2 


-  1.0  + 


(A20) 
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TABLE  I.-  MODEL  GEOMETRY 


Fuselage : 

Moment  reference  center: 
x . 

y . 

z . 

Length  . 

Rotor: 

Hub  coordinates: 

x . . 

y . 

z . 

Number  of  blades  .  .  . 

7/  deg  . 


hh  . 

Root  cutout  . 

Chord,  m  . 

Radius,  R,  m  .... 

Twist,  deg  . 

.  2 

Flapping  inertia,  kg-m 

Solidity  . 

Airfoil  section  .  .  . 


.  0 . 690R 

.  OR 

.....  OR 

.  2. OR 

.  0.69CR 

.  OR 

.  0.274R 

.  4 

.  0.009R 

.  0.034R 

.  0.090R 

.  0.20R 

.  0.108 

.  1.574 

.  -8.0 

.  0.653 

.....  0 • 0871 

NASA  RC-10-(B)  MOO 2 
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Figure  1.-  Axes  and  sign  conventions 


View  from  rear 

(b)  Rotor  system. 
Figure  1.-  Concluded. 


Figure  4.-  Wind-tunnel  model  installed  in  test  section  of  Langley  V/STOL  Tunnel. 


Figure  6.-  Isolated  rotor-induced  velocities 


Figure  9.-  Lateral  pressure  distributions  at  four  fuselage  stations 


Left  side 


Left  side 


Figure  9.-  Continued 


Left  side 


Left  side 


left  side  Right  side 


Figure  9.-  Continued. 


Separated  regi 


Unseparated  Separated 

Figure  10.-  Separation  and  swirl  effects  on  fuselage  cross-section 

pressure  distribution. 
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